Probing Our Built-in Calculator: A Systematic Narrative Review of Noninvasive Brain Stimulation Studies on Arithmetic Operation-Related Brain Areas

This systematic review presented a comprehensive survey of studies that applied transcranial magnetic stimulation and transcranial electrical stimulation to parietal and nonparietal areas to examine the neural basis of symbolic arithmetic processing. All findings were compiled with regard to the three assumptions of the triple-code model (TCM) of number processing. Thirty-seven eligible manuscripts were identified for review (33 with healthy participants and 4 with patients). Their results are broadly consistent with the first assumption of the TCM that intraparietal sulcus both hold a magnitude code and engage in operations requiring numerical manipulations such as subtraction. However, largely heterogeneous results conflicted with the second assumption of the TCM that the left angular gyrus subserves arithmetic fact retrieval, such as the retrieval of rote-learned multiplication results. Support is also limited for the third assumption of the TCM, namely, that the posterior superior parietal lobule engages in spatial operations on the mental number line. Furthermore, results from the stimulation of brain areas outside of those postulated by the TCM show that the bilateral supramarginal gyrus is involved in online calculation and retrieval, the left temporal cortex in retrieval, and the bilateral dorsolateral prefrontal cortex and cerebellum in online calculation of cognitively demanding arithmetic problems. The overall results indicate that multiple cortical areas subserve arithmetic skills.


Introduction
Dehaene's triple-code model (TCM) of numerical cognition draws a theoretical picture of number processing, including mental arithmetic, in the human brain (Dehaene et al., 2003).
Based on different streams of research, such as experimental psychology, neuroimaging, and neuropsychological testing in humans, as well as single-unit recordings in monkeys, the TCM suggests three distinct codes for numbers represented in segregated neural populations in the parietal cortex.

Assumptions of the TCM
The first TCM assumption is the presence of quantity representation or magnitude code at the bilateral intraparietal sulcus' horizontal segment (hIPS) for understanding and manipulating numerical quantities (Dehaene et al., 2003).Indeed, patients with IPS lesions have deficits in subtraction, an operation thought to require actual manipulation (online calculation) of numerical quantities but preserved knowledge of rote-learned arithmetic facts (Dehaene and Cohen, 1997;Cohen et al., 2000;Dehaene et al., 2003).The same holds for functional magnetic resonance imaging (fMRI) studies in healthy subjects, subtraction (with some exceptions of problems that can be stored in verbal memory such as 3 − 3 = 0), and two-digit addition and multiplication elicits left (Dehaene and Cohen, 1997;Kazui et al., 2000;Zago et al., 2001;Simon et al., 2002;Delazer et al., 2003Delazer et al., , 2005;;Kawashima et al., 2004;Ischebeck et al., 2006;Faye et al., 2019), right (Cohen et al., 2000;Kazui et al., 2000;Prado et al., 2011), and bilateral IPS activation (Lee, 2000;Menon et al., 2000;Kong et al., 2005;Zago et al., 2008;Andres et al., 2011;Klein et al., 2013a,b).The right hIPS is thought to be involved in forming spatial/structural diagrams that aid in spatial manipulation of numbers to obtain the solution (Miller et al., 2012).
The TCM also suggests a verbal code or verbal representation of numbers at the left angular gyrus (AG), where numerals are coded as words or sequences of words with lexical, phonological, and syntactical representations.The model predicted that simple or single-digit problems, such as multiplication, could be solved using operands (e.g., 2 × 5) transcoded into verbal code (two times five), which would elicit the rote memory (two times five equals ten) of this operation (Arsalidou and Taylor, 2011).Accordingly, left AG lesions impair simple addition, multiplication, and division but spare subtraction skills (Warrington, 1982;Lampl et al., 1994;Dehaene and Cohen, 1997;Cohen et al., 2000;Lee, 2000).In some cases, impaired multiplication with spared subtraction co-occurs with aphasia, indicating the dependency of multiplication on a language-based system of number processing (Dehaene and Cohen, 1997;Cohen et al., 2000;Dehaene et al., 2003).Multiplication engages the left AG because it is learned by rote memorization of multiplication tables in schools and, therefore, has strong associations between problems and solutions that are stored in long-term memory as arithmetic facts (Dehaene et al., 2003;Wassermann et al., 2008;Klein et al., 2013a,b;Faye et al., 2019).Single-digit additions (or subtraction) with a sum below 10 (e.g., 2 + 2 = 4) are also stored and retrieved from memory using the verbal code (but see for a different view).Increased left AG BOLD activity is observed during long-term memory retrieval (Menon et al., 2000;Stanescu-Cosson et al., 2000;Grabner et al., 2009;Jost et al., 2009;Prado et al., 2011;Klein et al., 2013a,b;De Visscher et al., 2015;Soylu and Newman, 2015) or when quantity-based processes shifted to automatic arithmetic fact retrieval (Ischebeck et al., 2006;Bloechle et al., 2016).Incorrect or "confusion" equations in which the proposed answer was true for a related operation (e.g., 9 × 6 = 15) also increased the left AG activation because they activate related arithmetic facts stored in verbal long-term memory (Grabner et al., 2013).The right AG activation was also reported and thought to be due to visuospatial attending processes crucial during fact retrieval (Grabner et al., 2009;Arsalidou and Taylor, 2011).
The model's third assumption is the presence of an Arabic code or visual representation of numbers at the bilateral posterior superior parietal lobule (PSPL).This code is thought to be used in accessing the other two systems through the visual modality and monitoring complex arithmetical calculations, such as those that require borrowing and carrying procedures, and magnitude code-independent operations, such as parity judgments (Dehaene et al., 1993;Wassermann et al., 2008).Indeed, there is significant bilateral PSPL activation during two-digit subtraction task performance or when two operations are performed instead of one (Lee, 2000;Menon et al., 2000;Andres et al., 2011).PSPL activation is also found to be problem size-, strategy-, and operands number-dependent (Menon et al., 2000;Stanescu-Cosson et al., 2000;Delazer et al., 2005;Zhou et al., 2007;Andres et al., 2011).Nonetheless, the PSPL's function is considered not specific for numerical computations since manual, visuospatial (e.g., saccade tasks), grasping, pointing, or linguistic tasks also activate it (Simon et al., 2002;Knops et al., 2009).
In recent years, theories that deviate from the original conception of TCM have been proposed (e.g., symbol-to-referent mapping hypothesis (Walsh, 2003), and updates of the model were suggested.This movement has been partly influenced by the "number neurons" discovery in the IPS and prefrontal cortex (PFC) of nonhuman primates (Nieder, 2016).The findings in animals are mirrored in human brain imaging studies; numerical processing robustly activates the inferior and superior parietal lobules, IPS, inferior and middle frontal gyrus, and cingulate gyrus (Arsalidou and Taylor, 2011;Nieder, 2016;Hawes et al., 2019).The nonparietal areas are now considered crucial for arithmetic proficiency as they serve domain-general processes, including working memory, episodic and semantic memory, and executive control (Garcia-Sanz et al., 2022).

Noninvasive brain stimulation (NIBS)
NIBS techniques, such as transcranial magnetic stimulation (TMS) and transcranial electrical stimulation (tES) that can modulate neuronal excitability and induce neuroplastic effects, have rekindled interest in exploring the causal role of the parietal cortex for arithmetic processing.Traditionally, these methods were classified as either "facilitatory" or "inhibitory" based on their effect on motor-evoked potentials (MEPs) elicited by the singlepulse TMS stimulation of the contralateral motor cortex (Boroojerdi et al., 2001;Paulus et al., 2008;Ziemann et al., 2008;Klomjai et al., 2015).The MEP amplitude depends on motor pathway neuron depolarization with the electrical current induced by the strong time-varying magnetic field; hence, it is considered an index of corticospinal excitability.Single-pulse TMS has a short-lived effect; however, when applied repeatedly (repetitive TMS or rTMS) at low frequencies (<1 Hz) or high frequencies (>5 Hz), it can induce a build-up of sustained corticospinal excitability suppression and facilitation, respectively (Pascual-Leone et al., 1994;Chen et al., 1997;Wassermann et al., 1998;Houdayer et al., 2008).Similar effects are produced with repeated application of bursts containing three TMS pulses of high frequency (50 Hz) repeated within the theta frequency range (Huang et al., 2005).These patterned rTMS protocols are known as theta burst stimulation (TBS) and reduce cortical excitability when applied continuously for 20 or 40 s (continuous TBS) and increase cortical excitability when 2 s trains are repeated every 10 for 190 s (intermittent TBS).The aftereffects of rTMS mimic synaptic plasticity in animal models, specifically, namely, long-term depression (LTD) and long-term potentiation (LTP; Urban et al., 2002;Gersner et al., 2011).For instance, the long-lasting potentiation and depression of corticospinal excitability induced by TBS can be blocked by a N-methyl-D-aspartate (NMDA) receptor antagonist (Wankerl et al., 2010), as well as the long-lasting potentiation elicited by high-frequency rTMS (Schwenkreis et al., 2005).Meanwhile, low-frequency rTMS enhances the amplitude and duration of γ-aminobutyric acid (GABA-B)-mediated inhibitory postsynaptic potentials (Casula et al., 2014).
For tES, low-intensity (1-2 mA) electrical current is applied to the head using surface electrodes to modulate cortical neural activity (Reed and Cohen Kadosh, 2018).The most widely used method, transcranial direct current stimulation (tDCS), involves the continuous application of direct electrical currents, typically for 10 to 30 min.The concurrent effects of the bipolar tDCS stimulation of the motor cortex are akin to a "somatic doctrine" of membrane polarization.The inward current flow (injection of cations) from the anode electrode to the brain leads to the hyperpolarization of apical dendrites and then to the depolarization (excite) of the neuronal soma and axon hillock, whereas the outward current flow from the brain to the cathode causes the depolarization of apical dendrites (accumulation of negative charges on the outer surface of the membrane) and then the hyperpolarization (inhibition) of the neuronal soma (Nitsche et al., 2008;Stagg and Nitsche, 2011;Lefaucheur and Wendling, 2019).Meanwhile, the polarity-dependent aftereffects of tDCS on motor cortical excitability depend on a similar synaptic mechanism demonstrated in animals for glutamatergic LTP and LTD (Bear and Malenka, 1994).The anodal tDCS elicited potentiation, and the cathodal tDCS-elicited depression of MEP amplitudes beyond stimulation was shown to be NMDA-and calcium channel-dependent (Liebetanz et al., 2002;Nitsche et al., 2003Nitsche et al., , 2004)).
There are also tES methods involving the application of rhythmic alternating current at a single frequency (transcranial alternating current stimulation or tACS) or at multiple changing frequencies ranging from 0.1 to 640 Hz and intensity (transcranial random noise stimulation or tRNS) or between 1 and 5 Hz (pulsed low-amplitude alternating electrical current or tPCS).For tACS, the online effect is attributed to the entrainment of endogenous oscillatory activity by the exogenously applied oscillating sinusoidal current in a frequency-specific (phase-to-phase locking) or cross-frequency (phase-amplitude coupling) manner (Antal and Paulus, 2013;Riddle et al., 2021).Entraining implies that the phase and/or frequency of the brain oscillation is modulated to follow external stimulation (Herrmann and Strüber, 2017).For instance, tACS applied at a fixed frequency (10 Hz) has no robust effect on corticospinal excitability (Antal et al., 2008) but increases it when applied at individual participants' alpha peak frequencies determined with electroencephalography recordings (Fresnoza et al., 2018).Meanwhile, the evidence indicates that plasticity mechanisms are sufficient to explain the aftereffects of tACS (Vossen et al., 2015) since an NMDA receptor antagonist suppresses the effect of 20 Hz tACS on corticospinal excitability and beta oscillations (Wischnewski et al., 2019).
In contrast to tACS, tRNS is applied at a wide frequency spectrum, between 0.1 and 100 Hz (low-frequency tRNS), between 101 and 640 Hz (high-frequency tRNS), and between 0.1 and 640 Hz (full-spectrum tRNS).For corticospinal excitability, lowfrequency tRNS has no measurable effects, whereas high-frequency and full-spectrum tRNS increases it beyond the stimulation period (Terney et al., 2008;Moliadze et al., 2014).Studies have shown that the sodium channel blocker carbamazepine tends to inhibit corticospinal excitability after tRNS stimulation in humans (Chaieb et al., 2015).This supports the proposed mechanism behind tRNS online effects, which is the repetitive opening of the sodium channels, shown during the application of AC stimulation to rat hippocampal slices (Schoen and Fromherz, 2008).For the aftereffects, however, a plasticity mechanism is unlikely responsible because the partial NMDA receptor agonist D-cycloserine and the NMDA receptor antagonist dextromethorphan had no significant impact on the excitability increases after tRNS (Chaieb et al., 2015).Instead, the alternative proposed mechanism is stochastic resonance, where the introduction of white noise by tRNS boosts the synchronization of neural firing through the amplification of subthreshold oscillatory activity, which in turn reduces the amount of endogenous noise (Van Doren et al., 2014;Antal and Herrmann, 2016;Jaušovec and Pahor, 2017).Meanwhile, in tPCS, where the current is delivered through bilateral earlobe electrodes, the mechanism of action involves the modulation of the temporal cortex neuronal activity and frontotemporal functional connectivity by increasing interhemispheric coherence at low-range frequencies (Datta et al., 2013;Morales-Quezada et al., 2015).
In the study of higher-order cognitive functions such as arithmetic, NIBS can be applied "online" or simultaneously with task execution to modulate task-relevant ongoing neuronal activity.The injected current can serve as additional neural noise that can be beneficial if it is too weak to overwhelm the neural network and could actually boost the overall background activity the signal of interest is riding on.In some cases, however, the noise is so robust that it briefly interferes with the brain's ability to selectively disrupt the processing of distracting stimulus elements (Luber and Lisanby, 2014).In the case of tACS or rTMS, the stimulation could directly change neural dynamics by resetting, driving, or enhancing rhythmic neuronal firing that has functional relevance (Luber and Lisanby, 2014).For example, in the IPS, which is a part of the dorsal attention network involved in endogenous shifts of attention, individual alpha peak frequency-tuned tACS caused an attentional bias in a spatial cueing task (Kemmerer et al., 2022).This finding not only supports the functional role of alpha oscillations in visuospatial attention but also suggests that spatial attentional processing, which is also employed in solving arithmetic problems, is a highly probable process disrupted or enhanced by NIBS during online stimulation (Masson and Pesenti, 2014).
NIBS is also applied "off-line," either before a task to induce up-or downregulation of intrinsic neuronal excitability that can affect task performance afterward, or in some cases, after a task to modulate consolidation processes in a learning paradigm (Hartwigsen and Silvanto, 2023).As the extension of effects beyond the stimulation period is attributed to neuroplasticity (except for tRNS), the potentiation and depression of synaptic strength would be the likely mechanism of arithmetic task performance improvement or impairment, respectively.An alternative mechanism put forward by Luber and Lisanby (2014) for the enhancing effect of TMS is called the "addition-by-subtraction" principle.This mechanism suggests that TMS might produce cognitive enhancement by disrupting or inhibiting an inessential or less essential but competing part of one or more functional brain networks involved in a task (Luber and Lisanby, 2014).For example, applying an inhibitory paradigm (e.g., 1 Hz rTMS) to brain areas contralateral to those relevant to the task might release interhemispheric inhibition.

The present review
In the last decades, NIBS studies have provided us with a considerable amount of data to affirm and even challenge the TCM assumptions.However, the TCM has not been fully tested against the findings of all NIBS-arithmetic studies available to date because previous reviews focused on different aspects of numerical cognition (e.g., numerosity and magnitude representation), only one method (e.g., tES or TMS), or only one assumption (Sarkar and Cohen Kadosh, 2016;Schroeder et al., 2017;Faye et al., 2019;van Bueren, 2021a;Garcia-Sanz et al., 2022;Lazzaro et al., 2022).Therefore, this review comprehensively synthesized results from all NIBS studies involving simple and complex symbolic arithmetic tasks in healthy adults and patient groups to test the TCM assumptions.

Materials and Methods
This review was conducted according to the accepted guidelines for the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (Fig. 1; Page et al., 2021).
Eligibility criteria and quality assessment.For the selection of relevant articles, the study had to (1) include a sample of healthy subjects or patients solving arithmetic problems before and after, or during TMS, tDCS, tACS, tRNS, and tPCS stimulation of parietal and nonparietal brain areas; (2) be peer-reviewed; (3) be accessible (full-text) and written in English; (4) have a single-or double-blinded, randomized, sham-controlled, and crossover or parallel study design; (5) include participants ≥18 years of age; and (6) be carried out in accordance with the Helsinki Declaration.We excluded editorials, commentaries, meeting or conference abstracts, animal studies, and NIBS studies concerning nonarithmetical aspects of numerical cognition, which were discussed extensively elsewhere (Dehaene et al., 2003;Berch and Mazzocco, 2007;Wassermann et al., 2008;Looi and Cohen Kadosh, 2016;Sarkar and Cohen Kadosh, 2016;Schroeder et al., 2017;Faye et al., 2019;Hawes et al., 2019;van Bueren, 2021b).Eligible full-text articles were screened for study design, blinding method, number of participants and mean age, stimulator model, stimulation parameters, target cortical areas, localization methods, arithmetic and control tasks, task time (online vs off-line), and the main findings.

Results
The initial search terms identified 103 potentially relevant articles published over the last 22 years .The screening process removed 63 articles, mainly NIBS studies concerning nonsymbolic arithmetic and magnitude tasks, as well as duplicate copies and studies with a nonretrievable full manuscript.One study in children (Looi et al., 2017) was excluded because the neurodevelopmental perspectives of arithmetic skills are not considered in the TCM, and the manner in which a child's brain engages in arithmetic might differ from that of an adult.For instance, greater PFC activation has been observed in younger children during arithmetic problem-solving, indicating greater cognitive and working memory demands during the early stages of skill acquisition (Rivera et al., 2005;Menon and Chang, 2021).In total, 37 original studies conducted with 1,219 participants met the inclusion criteria and were included in the review.Tables 1-3 provide an overview of each study's stimulation parameters.Thirty-three studies were conducted on healthy adults (n = 1,157) and four on patient groups (n = 62).Across the studies, the targets were the IPS (19), the AG (7), the PSPL (2), supramarginal gyrus (SMG), and nonparietal areas (24).The studies used various methods: tDCS (21), followed by TMS (8), tRNS (6), tACS (2), and tPCS (1).The primary outcome measures are reaction time (RT) and accuracy or error rates.None of the reviewed studies reported stimulation adverse or side effects.The following section describes the studies' main findings.

Effects of IPS stimulation on arithmetic performance
The TCM proposes the engagement of bilateral IPS whenever a task requires access to number magnitude representation and calculation procedure (Dehaene et al., 2003).This assumption seems to be valid since RTs for subtraction problems are prolonged by right or left hIPS 10 Hz rTMS (Andres et al., 2011;Montefinese et al., 2017) and for difficult addition problems by single-pulse TMS (Salillas et al., 2012).10 Hz rTMS and single-pulse TMS of the right or left ventral portion of the IPS (vIPS) also prolonged RTs for addition and subtraction and difficult multiplication problems, respectively (Salillas et al., 2012;Montefinese et al., 2017).Notable is the hemispheric interference asymmetry.In the Montefinese study, more interference was found for right than left hIPS stimulation for addition and subtraction and stronger for right hIPS than right vIPS stimulation for addition.In the Salillas study, interference for multiplication, which is related to problem size, was also stronger for the right vIPS.Left hIPS 1 Hz rTMS did not impair subtraction task performance, also indicating the right hemispheric dominance for online calculation (Fresnoza et al., 2020).Results from tES studies are also broadly consistent with the first assumption.For example, participants who received 1 and 2 mA left IPS anodal tDCS have shorter calculation times and higher scores than the control group (Houser et al., 2014(Houser et al., , 2015)).Left but not bilateral IPS anodal tDCS also shortened RTs for a two-digit subtraction task and increased accuracy in a magnitude comparison task (Hauser et al., 2013).Meanwhile, bilateral IPS tDCS with left anode-right cathode montage (but not the reverse montage) shortened RTs for two-digit times one-digit multiplication problems (e.g., 64 × 7) in participants with left calculation dominant hemisphere (Kasahara et al., 2013).In a two-digit addition verification task, simultaneous bilateral IPS anodal or cathodal stimulation also reduced and increased the distractor distance effect, respectively (Klein et al., 2013a,b), while right IPS anodal but not cathodal stimulation prolonged latency increases for carry operations on the same task (Artemenko et al., 2015).
The results of tRNS studies, on the other hand, point to a beneficial effect of IPS stimulation for enhancing arithmetic proficiency through learning rather than facilitating actual numerical task performance.As proof of this, bilateral IPS high-frequency tRNS shortened RTs and improved accuracy in two-digit subtraction verification tasks during testing and retesting regardless of problem novelty (Pasqualotto, 2016).Coupling training with bilateral IPS full-spectrum tRNS also stabilized the training-induced performance improvement, although the transfer effect was somewhat limited to related tasks or modalities (Cappelletti et al., 2013(Cappelletti et al., , 2015;;Popescu et al., 2016).In the Popescu study, 5 d training with drill and calculation problems coupled with bilateral PFC high-frequency tRNS (first 3 d) and IPS high-frequency tRNS (last 2 d) consequently improved RT's latency and distribution for calculation problems in the difficult condition during training and accuracy in the easy and new conditions during testing.In the Cappelletti studies, 5 d training of numerosity discrimination concurrent with bilateral IPS full-spectrum tRNS led to a longterm (16 weeks) improvement in the training task across age groups.However, the training had an age-dependent effect on the transfer effect to untrained related tasks: space and time discrimination was improved and impaired in young and old participants, respectively.Meanwhile, training did not significantly affect untrained and unrelated tasks, including simple and complex symbolic arithmetic tasks.
There are also results at odds with the first assumption of the TCM.Some studies showed an impairment of arithmetic fact retrieval under IPS stimulation: longer RTs were observed in one-digit multiplication during left but not right hIPS single-pulse TMS (Salillas et al., 2012) and in two-digit addition during 10 Hz rTMS applied between the left AG and the posterior part of the IPS (Göbel et al., 2006).Meanwhile, in another study, rote-learned one-digit subtraction and multiplication were unaffected by left hIPS 1 Hz rTMS (Fresnoza et al., 2020).Another result that is at odds with the first assumption of the TCM is that left PPC theta tACS decreased calculation times in novel fact-learning problems, as the theta frequency band indexes fact retrieval from memory (Mosbacher et al., 2021).Other results in contrast to the first assumption are the absent effects of right IPS anodal tDCS on two-digit subtraction and magnitude comparison (Hauser et al., 2013(Hauser et al., , 2016) ) and the absent effects of left IPS anodal tDCS on one-digit subtraction and one and two-digit additions, as well as two-back task performance (Mosbacher et al., 2020), and of left PPC alpha tACS and anodal tDCS on the training and acquisition of a novel arithmetic procedure and novel arithmetic facts (Mosbacher et al., 2021).Bilateral IPS highfrequency tRNS also failed to modulate target identity, distractor distance, and carry-over effects in a two-digit addition verification task (Bieck et al., 2018).

Effects of AG stimulation on arithmetic performance
The second TCM assumption states that the left AG is involved in the verbally mediated retrieval of symbolic number knowledge and arithmetic facts (Dehaene et al., 2003;Wassermann et al., 2008;Sokolowski et al., 2017).This is corroborated by several results: 5 Hz rTMS of the left AG or at an area very close to it (left anterior superior temporal gyrus) evoked the highest error rate for single-digit multiplication (30%) and one-digit addition (35%), respectively (Maurer et al., 2015).Furthermore, left AG 1 Hz rTMS prolonged RTs for subtraction problems often solved by retrieval (e.g., 14 − 7), simple multiplication problems, and online calculation (Fresnoza et al., 2020).There was also no effect of right AG anodal tDCS on one-digit multiplication verification task performance (Clemens et al., 2013).However, several results indicate involvement of the left or bilateral AG in complex calculations: right AG 5 Hz rTMS increased the error rate for simple addition problems (Maurer et al., 2015), 10 Hz rTMS to either the right or left AG prolonged RTs for complex addition and subtraction problems (Montefinese et al., 2017), and 10 Hz rTMS between the left AG and IPS prolonged RTs for double-digit addition (Göbel et al., 2006).Although Göbel et al. attributed the effect to memory retrieval difficulties, this is unlikely because their participants reported mentally calculating the solutions.In addition, Grabner et al. (2015) showed that left AG cathodal tDCS coupled with arithmetic fact training prolonged RTs for complex multiplication and subtraction that persisted over 24 h but only for trained problems.In the same study, left AG anodal tDCS and training improved subtraction but not multiplication learning (Grabner et al., 2015).Left AG anodal tDCS also improved response latencies in large operand addition problems assumed to be solved by calculation (two-digit/two-digit carry problem with addends ranging from 12 to 29) but decreased solution rates for small operand addition problems assumed to be solved by memory retrieval (one-digit/one-digit problems with addends between 2 and 8; Rütsche et al., 2015).

Effects of PSPL stimulation on arithmetic skills
The third TCM assumption suggests a supportive role for the PSPL for the other two brain areas because it is assumed to hold the visual representation of numbers (Dehaene et al., 1993;Wassermann et al., 2008).However, the two rTMS studies that directly stimulated the bilateral PSPL (Table 1) failed to demonstrate a significant effect on single-digit addition, subtraction, multiplication, and division task performance (Andres et al., 2011;Maurer et al., 2015).

Effects of SMG and nonparietal area stimulation on arithmetic skills
Dehaene and colleagues also outlined domain-general processes interfacing with numerical cognition, such as executive function, spatial attention, and working memory (Arsalidou and Taylor, 2011;Schroeder et al., 2017;Siemann and Petermann, 2018).Against this background, several studies that targeted nonparietal and parietal areas outside those postulated by the TCM enhanced arithmetic proficiency.

SMG
The neuroimaging literature associated the right SMG with spatial working memory and attention shifting, particularly for cognitively demanding arithmetic tasks (Silk et al., 2010;Abd Hamid et al., 2011;Faye et al., 2019).Indeed, 10 Hz rTMS elicited a lateralized (right SMG > left SMG) detrimental effect on two-digit subtraction but not addition problems, which the authors attributed to the former being relatively less automated and more difficult to solve than the latter (Montefinese et al., 2017).Consistently, compared with all cortical targets, the highest error rate (11%) for single-digit subtraction was elicited at the right posterior SMG by 5 Hz rTMS (Maurer et al., 2015).Meanwhile, for the left SMG, some neuroimaging studies suggested a role in arithmetic fact retrieval (Faye et al., 2019).Ten Hz rTMS of the anterior part of the left SMG prolonged RTs in a consumer-like discount calculation task involving multidigit mental addition or subtraction problems (Klichowski and Kroliczak, 2020), while when applied between the left SMG and anterior IPS, it impaired two-digit addition task performance (Göbel et al., 2006).

PFC
Arithmetic tasks involving procedural strategies recruit the frontal lobe, particularly the dorsolateral region of the PFC (DLPFC; Grabner et al., 2009;Arsalidou and Taylor, 2011;Hawes et al., 2019;Mosbacher et al., 2020).Interestingly, it is also one of the core brain regions showing elevated activity in calculation prodigies (Pesenti et al., 2001;Krause et al., 2019).However, instead of a specialized function in arithmetic, the PFC is ascribed a role in more general functions such as top-down attention, executive functions, and working memory, the cognitive space where number representations can be stored and manipulated (e.g., choice of strategy and planning; Arsalidou and Taylor, 2011;Klein et al., 2013aKlein et al., ,b, 2016;;Matejko and Ansari, 2015;Schroeder et al., 2017;Siemann and Petermann, 2018;Hawes et al., 2019).
A consistent finding of the studies reviewed here was improved arithmetic proficiency by unilateral or bilateral tES stimulation of the dorsolateral region of the PFC, particularly for cognitively demanding tasks (Table 2).For instance, left PFC anodal tDCS increased accuracy, shortened RTs, and decreased RT variability on the Paced Auditory Serial Subtraction Task (PASST) but had no effect on the less challenging Paced Auditory Serial Addition Task (PASAT; Pope et al., 2015).Meanwhile, in the more difficult adaptive PASAT (A-PASAT), left PFC anodal tDCS improved calculation speed and reduced emotionally salient stimuli-elicited attentional biases (Plewnia et al., 2015), as well as improved calculation accuracy when coupled with a complex three-back but not one-back training task (Gill et al., 2015).Left PFC anodal tDCS also shortened RTs for two-digit subtraction problems but had no effect on single-digit addition and subtraction problems, as well as on easy two-back working memory tasks (Mosbacher et al., 2020).Similarly, high-definition left PFC cathodal tDCS improved complex matchstick arithmetic task performance, which the authors attributed to the successful inhibition of previously learned rules that constrain the participants' ability to change problem representations (Luft et al., 2017).
Combining sinusoidal current stimulation with training also benefitted arithmetic performance and learning.Bilateral PFC high-frequency tRNS and rote memorization (Days 1-3) and bilateral IPS high-frequency tRNS and calculation training (Days 4-5) shortened RTs for single-and double-digit addition, subtraction, and difficult multiplication trials (with few problem repetitions) compared with easy trials (with many problem repetitions) during testing (Popescu et al., 2016).Bilateral PFC high-frequency tRNS coupled with algorithmbased calculation and memory retrieval-based drill training for 5 consecutive days also shortened RTs in both tasks and elicited sustained calculation improvement 6 months after training (Snowball et al., 2013).Left PFC theta tACS and 25 min arithmetic operation training of procedural and fact knowledge also accelerated the decrease in calculation times for fact knowledge (Mosbacher et al., 2021).On the other hand, nontraining studies indicate that bilateral PFC high-frequency tRNS is more effective during actual numerical performance testing than learning, as shown by the increased distractor distance effect in the easy and difficult trials of a double-digit addition verification task (Bieck et al., 2018).Meanwhile, in a double-digit subtraction verification task, bilateral PFC or IPS highfrequency tRNS shortened RTs during stimulation (online) and increased accuracy 7 d after stimulation (off-line; Pasqualotto, 2016).
There are also reports of null and state-dependent effects of frontal lobe stimulation on arithmetic performance.Left PFC anodal tDCS or alpha tACS did not affect learning novel arithmetic procedures and facts (Mosbacher et al., 2021), while concurrent right MFG and left IPS individually tuned HD theta tACS did not improve one-digit addition and two-digit subtraction task performance (Zhang et al., 2022).Interestingly, bilateral PFC fullspectrum tRNS failed to enhance the performance of a calculation prodigy in a complex multiplication task and reduced task accuracy of postgraduate students with highly standardized mathematics scores in simplified multiplication problems, indicating ineffective or deleterious effects for individuals with already high-performing arithmetic functioning (Krause et al., 2019).Correspondingly, anode-left and cathode-right PFC tDCS prolonged and shortened RTs for simple arithmetic decision tasks in low-and high-math anxiety individuals, respectively (A Sarkar et al., 2014).

Temporal cortex
The temporal cortex, including the AG, is part of the language-processing regions of the left hemisphere and serves as an additional storage of semantic knowledge (Andres et al., 2011;Siemann and Petermann, 2018).The AG is crucial for verbally mediated arithmetic fact retrieval, as demonstrated in a patient with left parietotemporal hemorrhage showing selective acalculia for addition, multiplication, and division but an intact ability to subtract (Lampl et al., 1994), and in a patient with a left superior temporal cortex lesion exhibiting a severe selective impairment in multiplication and division tasks (van Harskamp et al., 2005).In healthy participants, virtual lesions induced by 5 Hz rTMS elicited a higher error rate when applied on the left anterior STG (35%) than on the right posterior STG (8%) for simple addition tasks (Maurer et al., 2015).Meanwhile, bilateral earlobe tPCS did not affect accuracy for simple (e.g., 7 − 4) and difficult (e.g., 844 − 385) subtraction problems, ruling out the possibility that the temporal cortex performs online calculation procedures (Morales-Quezada et al., 2015).

Cerebellum
The cerebellum is crucial in developing early mathematical ability because of its role in sequence or pattern detection and optimization through constant error correction (Vandervert, 2017).In arithmetic performance in adults, the cerebellum is associated with various verbal working memory processes, including automated number manipulation, as reflected by the increase in BOLD signal at the cerebellar cortex and prefrontal and parietal cortices during PASAT task performance (Hayter et al., 2007).In the only NIBS study that stimulated the cerebellum, right cerebellar cathodal tDCS increased accuracy and decreased mean verbal RT and variability in the PASST task but not the PASAT task (Pope and Miall, 2012).The authors argue that these results are consistent with the notion that inhibiting the cerebellum or disinhibiting the PFC can help release cognitive resources needed for more difficult tasks (Pope et al., 2015).

NIBS studies in patient groups
In clinical populations, the results of few arithmetic studies were somewhat inconsistent (Table 3).In a brain tumor patient, 5 Hz rTMS neuronavigation (nrTMS)-based preoperative cortical mapping of simple arithmetic processing elicited the highest error ratio for onedigit addition at the left posterosuperior frontal gyrus and left AG, partially supporting the second TCM assumption.In contrast, the highest error ratio for multiplication, subtraction, and division was elicited in frontal lobe areas, indicating reliance on working memory processes.Importantly, the resection of nrTMS arithmetic-positive cortical sites disrupted arithmetic skills as measured by the standardized number-processing and calculation task (NPCT; Ille et al., 2018).In the case of two individuals with developmental dyscalculia, left anode-right cathode bilateral PPC tDCS coupled with numerical task training (Day 1, artificial digit task with Gibson figures; Days 2-6, artificial digit task, numerical Stroop and number line task) led to numerical Stroop and number line task performance improvements during retesting (Iuculano and Cohen Kadosh, 2014).Although arithmetic proficiency was not measured directly here, this study clarifies the PPC's and probably the IPS's function in manipulating numerical representations crucial for solving complex arithmetic problems.
Two other studies reported null results.In depression patients, although 2 week A-PASAT training (five sessions per week) with bilateral PFC tDCS (left anode-right cathode) reduced maladaptive ruminative thoughts, it failed to affect depression severity and arithmetic performance (Vanderhasselt et al., 2015).Meanwhile, in a patient with cerebral hypoxic ischemia-induced bilateral parietal lesions and Gerstmann syndrome, left PFC anodal tDCS immediately improved phonemic and semantic fluency (Sreeraj et al., 2019).Six months after hospital discharge, a 10 d consecutive left PPC anodal tDCS improved several cognitive task performances (Table 3), but not one-digit multiplication and division and all three-digit arithmetic task performances (Sreeraj et al., 2019).

Discussion
The present review was compiled to test the TMC assumptions against the results of 37 arithmetic studies that stimulated the IPS, AG, and PSPL, as well as other cortical areas with TMS and tES.Here, we will fit the overall results with each assumption, determine whether they support them, address results inconsistencies, and give a conclusion.
Concerning the second TCM assumption, the reviewed studies provide mixed results.Single-digit addition and multiplication performance disruption by left AG rTMS (Maurer et al., 2015;Fresnoza et al., 2020) and tDCS (Clemens et al., 2013) provided evidence in favor of the second TCM assumption.In contrast, the disruption of complex arithmetic performance does not support the assumed role of the left AG in arithmetic fact retrieval (Göbel et al., 2006;Grabner et al., 2015;Rütsche et al., 2015;Fresnoza et al., 2020).It could be argued that complex arithmetic problems are solvable by retrieval, namely, by breaking them down into simpler, retrievable problems (Dehaene and Cohen, 1995;Fayol and Thevenot, 2012;Siemann and Petermann, 2018).However, some new theories offer further explanations.The symbol-to-referent mapping hypothesis suggests that the left AG can match quantities to Arabic digits (Price and Ansari, 2011;Bloechle et al., 2016), while the attention to memory model proposes that the left AG is a flexible interface that adjusts and adapts attentional demands: top-down for magnitude manipulations and bottom-up for fact retrieval (Cabeza et al., 2008(Cabeza et al., , 2012;;Klein et al., 2013a,b;Bloechle et al., 2016).In other words, the AG could subserve fact retrieval for rote-learned arithmetic problems and online calculation for complex problems.The only drawback is that participant's trial-by-trial strategies are typically not assessed to confirm such an assumption.
For the third assumption, namely, that the PSPL, bilaterally, plays a role in arithmetic, two studies showed the absence of bilateral PSPL rTMS effect on single-digit addition, subtraction, multiplication, and division task performance (Andres et al., 2011;Maurer et al., 2015).This indicated that the PSPLs have no obligatory role in arithmetic processing and may only aid numeric computations by serving as a visuospatial medium for estimations (Chochon et al., 1999;Nieder and Dehaene, 2009) or by orienting and shifting attention toward the left or right side of a mental numerical continuum, helpful for the intuition of removing or adding quantities (Dehaene et al., 2003;Knops et al., 2009;Andres et al., 2011).So far, there is insufficient evidence for a definitive answer and that the third assumption requires further investigation.

Other cortical targets
The close examination of studies involving brain areas excluded from the main TCM assumptions showed that they also contribute to arithmetic processing.Results suggest the involvement of the right SMG for problems requiring online calculation (Maurer et al., 2015;Montefinese et al., 2017) and the left SMG for problems that require retrieval of answers from long-term memory (Göbel et al., 2006;Klichowski and Kroliczak, 2020).These result patterns are consistent with the meta-analysis of neuroimaging studies showing a clear link between the left SMG and arithmetic facts, as well as an association between the right SMG and approximate nonsymbolic magnitude processing (Faye et al., 2019).Interestingly, the interhemispheric differences in the SMG function were similar to those of the AG (Montefinese et al., 2017).It could be that the anatomical proximity of the SMG and AG is responsible for these results.
The evidence from frontal lobe stimulation points to three key findings.First, PFC stimulation reduces the impact of negative emotions, anxiety, distractors, and previously learned rules on arithmetic performance (Sarkar et al., 2014;Plewnia et al., 2015;Luft et al., 2017;Bieck et al., 2018).The consensus is that PFC stimulation reduces interference from these impediments, preserves working memory resources, and reduces attentional biases.Second, the stimulation has different effects on arithmetic performance for different groups of individuals.Bilateral PFC full-spectrum tRNS impairs individuals with high-performing arithmetic skills (Krause et al., 2019), while bilateral PFC tDCS improves and impairs the performance of individuals with high or low mathematics anxiety, respectively (Sarkar et al., 2014).
Third, the PFC, mainly the dorsolateral region, supports complex arithmetic task performance.It has been shown that left PFC tDCS (Gill et al., 2015;Plewnia et al., 2015;Pope et al., 2015;Luft et al., 2017;Mosbacher et al., 2020), as well as bilateral PFC highfrequency tRNS with training (Snowball et al., 2013;Popescu et al., 2016) or without training (Pasqualotto, 2016;Bieck et al., 2018), had beneficial effects on complex arithmetic task performance.These results are consistent with the idea that complex calculations have a higher demand for domain-general processes, including working memory (e.g., to maintain intermediate results) and attention, as they involve numerical decisionmaking and require many steps to generate a correct answer (Arsalidou and Taylor, 2011;Siemann and Petermann, 2018;Klichowski and Kroliczak, 2020).A hierarchical organization of the frontal gyri in arithmetic was also proposed.The inferior frontal gyrus is assumed to be involved in processing simple numerical tasks with few storage or procedural requirements, the middle frontal gyri in managing several cognitive procedural steps (e.g., carrying a number in two-digit addition), and the medial and superior frontal gyri in generating strategies to solve complex multistep problems (Arsalidou and Taylor, 2011).
The role of the temporal cortex and cerebellum in arithmetic warranted further explorations.The evidence from a few studies reported here suggests a role for the left temporal cortex in arithmetic fact retrieval since 5 Hz rTMS of the left anterior STG induced a high error rate for a simple addition (Maurer et al., 2015) and that tPCS did not affect simple and difficult subtraction problems that require online calculation procedures (Morales-Quezada et al., 2015).Meanwhile, the only study that stimulated the right cerebellum with cathodal tDCS revealed its supportive role in domain-general processes such as working memory, particularly for cognitively demanding arithmetic tasks such as PASST (Pope and Miall, 2012).

Limitations of the literature
The reviewed articles have potential limitations that must be mentioned.First, suitable control conditions are missing in some studies.This means that the domain specificity of the results for arithmetic is unclear.Second, divergent results can be linked to task heterogeneity.Subtraction and multiplication are commonly used to tap into the brain networks subserving the mental manipulation of numerical quantities and arithmetic fact retrieval, respectively.This is problematic since arithmetic problems can be solved using both strategies.However, except for Fresnoza et al. (2020) and Mosbacher et al. (2020), no other study considered the impact of trial-by-trial solution strategies on arithmetic performance.Third, the modulatory effects of NIBS on cognitive function are sensitive to a range of methodological-and subject-specific factors that are rarely considered in many studies [for an extensive review, see Ridding and Ziemann (2010), Paulus et al. (2016), Sarkar and Cohen Kadosh (2016), Silvanto and Cattaneo (2017), and Hartwigsen and Silvanto (2023)].For instance, stimulation intensities are typically set to values known to inhibit and facilitate corticospinal excitability.However, studies have shown that distinct TMS intensity ranges can induce distinct effects: low intensities facilitate neural activity and behavior, whereas high intensities induce suppression (Silvanto and Cattaneo, 2017).Furthermore, these ranges are shifted by changes in neural excitability.For instance, a TMS intensity, which normally induces suppression, can have a facilitatory effect if the stimulated neurons are being inhibited by ongoing task-related processes or preconditioning (Silvanto and Cattaneo, 2017).Therefore, stimulation intensity and brain-state interactions could explain why NIBS methods' distinction as either inhibitory or facilitatory does not always map exactly onto the respective behavioral outcomes (Hartwigsen and Silvanto, 2023).
The notion that NIBS does not work for everyone can also be accounted for by interindividual variability in the response to stimulation protocols.For instance, anodal tDCS only increased corticospinal excitability in individuals who are sensitive to single-pulse TMS (Labruna et al., 2016).Similarly, anodal tDCS improved response inhibition (measured by a stop signal task) in individuals with high brain connectivity but not in those with low connectivity (Li et al., 2019).On the other hand, the impact of left PFC cathodal tDCS was predicted by the baseline GABA-glutamate ratio.Individuals with higher level of inhibition (more GABA relative to glutamate) experienced more disruption in response selection paradigms (Filmer et al., 2019).It is interesting to note that the effect of anodal tDCS is also influenced by attention allocation.Participants who were engaged in a cognitive task during motor cortex stimulation did not experience an increase in corticospinal excitability (Antal et al., 2007).This highlights the need to identify and control factors that may affect the response to NIBS.While controlling the neuronal state is difficult, priming measures and ensuring participants are attentive yet relaxed can help reduce outcome variations.However, factors such as age, gender, and genetic polymorphisms cannot be changed, so controlling them through careful experimental design is crucial (Huang et al., 2017).
Lastly, the lack of tES focality is another limitation in terms of site-specific targeting.TES effects are limited not only to cortical areas beneath the electrodes but also to areas affected by the electric field generated between them, even reaching deep regions (Bikson et al., 2010;Garcia-Sanz et al., 2022).Out of all the analyzed tES studies, only seven conducted electric field distribution simulations (Klein et al., 2013a,b;Artemenko et al., 2015;Grabner et al., 2015;Hauser et al., 2016;Luft et al., 2017;Bieck et al., 2018;Zhang et al., 2022).This means that in the remaining studies, the target areas may not be adequately reached and optimally stimulated, specifically during IPS stimulation, where the depth of current penetration is critical.Nevertheless, the studies that did use finite-element method-based electric field calculations showed activation peaks in the IPS higher than what is recommended (0.017 mA/cm 2 ) to modify cortical excitability by tDCS in humans (Klein et al., 2013a,b;Artemenko et al., 2015;Hauser et al., 2016;Bieck et al., 2018;Zhang et al., 2022).It is essential for future studies to perform subjectspecific electric field simulation to select the optimal current dose and pinpoint cortical areas that are likely to be activated by tES (Kashyap et al., 2021;Van Hoornweder et al., 2023).Techniques such as transcranial-focused ultrasound, transcranial temporal interference stimulation, and transcutaneous vagus nerve stimulation that offer superior spatial resolution and penetration depth compared with conventional stimulation paradigms could also potentially resolve the focality issue in future studies.Machine learning methods (e.g., personalized Bayesian optimization) can also be used to search, learn, and recommend neurostimulation parameters based on a participant's baseline arithmetic abilities (van Bueren, 2021b).Most importantly, further studies are warranted to better understand putative mechanisms underlying the effect of NIBS methods.This will help validate inferences on the causal role of a cortical area for arithmetic with absolute certainty.

Conclusions
This review provides a new perspective on TCM based on the results of arithmetic processing modulated by NIBS methods.Taken together, the results yielded three main conclusions.First, overall findings corroborated the first assumption of the TCM that the bilateral IPS is involved in numerical quantity manipulation.Second, the second assumption of the TCM that the left AG has an exclusive role in arithmetic fact retrieval cannot be confirmed because of mixed results.The third assumption about the role of the bilateral PSPL in the number domain could also not be supported due to the limited number of studies.Finally, collective results from studies that stimulated brain areas outside those postulated by the TCM showed that the bilateral SMG and bilateral DLPFC are also crucial for arithmetic performance.In contrast, although promising, the evidence from studies involving the temporal cortex and cerebellum is still insufficient.In summary, the available evidence from NIBS studies favors the existence of a wider cortical network rather than a one-to-one correspondence between numerical abilities and the specific parietal brain areas postulated by the TCM (Nieder, 2016).

Figure 1 .
Figure 1.Number of included studies following the PRISMA flow diagram for systematic reviews and metaanalysis.

Table 1 .
Summary of stimulation parameters and results of TMS studies in healthy participants -No control task • Left MFG: 45% ER for division • Right AG: 40% ER for subtraction • Left anterior STG: 35% for addition (8% on the right posterior STG) • Left AG: 30% ER for multiplication • Bilateral PSPL: no error recorder for all types of arithmetic tasks (Table continues.)

Table 2 .
Summary of stimulation parameters and results of TES studies in healthy participants • No effect on the control task(Table continues.)

Table 3 .
Table continues.)Summary of stimulation parameters and results of NIBS studies in patients